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Abstract. The constraint on the ADD model of extra dimensions coming from 
photon annihilation into Kaluza-Klein graviton in supernova cores is revisited. In 
the two photon process for a conservative choice of the core parameters, we obtain the 
bound on the fundamental Planck scale Af, ^1.6 TeV. The combined energy loss rate 
due to nucleon-nucleon brehmstrahlung and photon annihilation processes is rederived, 
which shows that the combined bounds add only second decimal place to M,. The 
present study can strengthen the results that are available in the current literature for 
the graviton emission from SN1987A which puts a very strong constraints on models 
with large extra dimensions for the case of n = 3 . 
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1. Introduction 

Stars are potential sources for weakly interacting particles such as neutrinos, gravitons, 
axions, and other new particles that can be produced by nuclear reactions or by thermal 
processes in the hot stellar interior. The solar neutrino flux is now routinely measured 
with such a precision that compelling evidence for neutrino oscillations has accumulated. 
The measured neutrino burst from supernova SN1987A has been used to derive many 
useful limits. Even when the particle flux can not be measured directly, the absence of 
visible decay products, notably x- or 7-rays, can provide important information. The 
properties of stars themselves would change if they lose too much energy into a new 
channel. This "energy-loss argument" has been widely used to constrain a long list of 
particles and there properties. All of this has been extensively reviewed [H [2] 

The extra dimensional scenario due to Arkani-Hamed, Dimopoulos and Dvali 
(ADD) [3] , model predicts a variety of novel signals which can be tested using table-top 
experiments, collider experiments, astrophysical or cosmological observations. It has 
been pointed out that one of the strongest bounds on models of extra dimensions comes 
from SN1987A [4]. Various authors have done calculations to place such constraints on 
the extra dimensions [5]- [in]. In this paper, we calculate the energy loss rate due to 
graviton emission from SN1987A by photon-photon annihilation and derive the bounds 
on extra dimensions. We combine the result with that of nucleon-nucleon brehmstralung 
process and derive the corresponding bound on large extra dimensions. 

Physically, there are two fundamental types of supernovae (SNe), based on what 
mechanism powers them: the thermonuclear SNe (Type I SNe) and the core-collapse 
ones (Type II SNe). The core-collapse SNe are the class of explosions which mark 
the evolutionary end of massive stars (M > 8Mq). Such stars have the usual onion 
structure with several burning shells, an expanded envelope, and a degenerate iron core 
that is essentially an iron white dwarf. The core mass grows by the nuclear burning at 
its edge until it reaches the Chandrasekhar limit. The collapse can not ignite nuclear 
fusion because iron is the most tightly bound nucleus. Therefore, the collapse continues 
until the equation of state stiffens by nucleon degeneracy pressure at about nuclear 
density (3 x lO^^'gcm"^). At this "bounce" a shock wave forms, moving outward and 
expelling the stellar mantle and envelope. The explosion is a reversed implosion, the 
energy derives from gravity, not from nuclear energy. Within the expanding nebula, a 
compact object remains in the form of a neutron star or perhaps sometimes a black 
hole. The kinetic energy of the explosion carries about 1% of the liberated gravitational 
binding energy of about 3 X 10^^ erg, 99% going into neutrinos. This powerful and 
detectable neutrino burst is the main astro-particle interest of core-collapse SNe. In 
core-collapse SNe only 10~^ of the total energy shows up as light, i.e. about 1% of 
the kinetic explosion energy, hence they are dimmer than SNe-Ia, and are not useful as 
standard candles. 

In the case of SN1987A, about 10^^ ergs of gravitational binding energy was released 
in few seconds and the neutrino fluxes were measured by Kamiokande [11] and 1MB [12] 
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collaborations. Numerical neutrino light curves can be compared with the SN1987A data 
where the measured energies are found to be "too low." For example, the numerical 
simulation in p!3] yields time-integrated values {E^J ^ 13MeV, {EpJ ^ 16MeV, and 
{E^^) ?a 23MeV. On the other hand, the data imply {EpJ = 7.5 Me V at Kamiokande 
and 11.1 MeV at 1MB [14J. Even the 95% confidence range for Kamiokande implies 
{EpJ < 12 MeV. Flavor oscillations would increase the expected energies and thus 
enhance the discrepancy [H]. It has remained unclear if these and other anomalies 
of the SN1987A neutrino signal should be blamed on small-number statistics, or point 
to a serious problem with the SN models or the detectors, or is there a new physics 
happening in SNe? 

Since we have these measurements already at our disposal, now if we propose some 
novel channel through which the core of the supernova can lose energy, the luminosity 
in this channel should be low enough to preserve the agreement of neutrino observations 
with theory. That is, 

^new channel ^ 10 CrgS S . (l) 

This idea was earlier used to put the strongest experimental upper bounds on the axion 
mass [15]. Here, we consider the emission of the higher- dimensional gravitons from the 
core. Once these particles are produced, they can escape into the extra dimensions, 
carrying energy away with them. The constraint on the luminosity of this process can 
be converted into a bound on the fundamental Planck scale of the theory, M*. The 
argument is very similar to that used to bound the axion-nucleon coupling strength 
[H [El [TTJ [18]. The 'standard model' of supernovae does an exceptionally good job of 
predicting the duration and shape of the neutrino pulse from SN1987A. Any mechanism 
which leads to significant energy-loss from the core of the supernova immediately after 
bounce will produce a very different neutrino-pulse shape, and so will destroy this 
agreement as demonstrated explicitly in the axion case by Burrows, Brinkmann, and 
Turner [18] . Raffelt has proposed a simple analytic criterion based on detailed supernova 
simulations [1]: if any energy-loss mechanism has an emissivity greater than 10^^ ergs 
g~ls~l then it will remove sufficient energy from the explosion to invalidate the current 
understanding of SNe II neutrino signal. 

2. Supernovae and constraints on large extra dimensions 

The most restrictive limits on come from SN1987A energy-loss argument. If large 
extra dimensions exist, the usual four dimensional graviton is complemented by a tower 
of Kaluza-Klein (KK) states, corresponding to new phase space in the bulk. The KK 
gravitons interact with the strength of ordinary gravitons and thus are not trapped in 
the SN core. During the first few seconds after collapse, the core contains neutrons, 
protons, electrons, neutrinos and thermal photons. There are a number of processes in 
which higher-dimensional gravitons can be produced. For the conditions that pertain 
in the core at this time (temperatures T ~ 30 — 70 MeV, densities p ~ (3 — 10) x 10^^ g 
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cm~^), the relevant processes are nucleon-nucleon brehmstrahlung, graviton production 
in photon fusion and electron-positron annihilation. 

In SNe, nucleon and photon abundances are comparable (actually nucleons are 
somewhat more abundant). In the following we present the bounds derived by various 
authors using nucleon-nucleon brehmstralung and in the next section we give a detailed 
calculation for photon-photon annihilation to KK graviton process. 

2.1. Nucleon- Nucleon brehmstralung 

This is the dominant process relevant for the SN1987A where the temperature is 
comparable to pion mass m^^ and so the strong interaction between nucleons is 
unsuppressed. This process can be represented as 

N + N^N + N + KK (2) 

where N can be a neutron or a proton and KK is a higher- dimensional graviton. 

The main uncertainty comes from the lack of precise knowledge of temperatures in 
the core: values quoted in the literature range from 30 to 70 MeV. For T = 30 MeV 
and p = 3 X 10^^ g cm~^, we list the results obtained by various authors. 
CuUen and Perelstein [5] 

n = 2, e = 6.79 x 10^^ x M^^ erg g-^ M, > 50TeV; (3) 

n = 3, e= 1.12 x 10^^ x M^^ erg g-^ M, > 4TeV. (4) 

Barger, Han, Kao and Zhang [6] 

n = 2, e = 6.7 X 10^^ x M^^ erg g-^ s-\ M, > 51TeV; (5) 

n = 3, e = 6.3 X 10^^ x M^^ erg g-^ s^\ M, > 3.6TeV. (6) 

Hanhart et. al. \T, "8] 

n = 2, e = 9.24 x 10^^ x M^^ erg g'^ s-\ M, > 31TeV; (7) 

n = 3, e = 1.57 X 10^^ x M^^ erg g-^ > 2.75TeV. (8) 

Hannestad and Raffelt P [ID] 

n = 2, e = 4.98 x 10^^ x M^^ erg g-^ s-\ M, > 84TeV; (9) 

n = 3, e = 1.68 x 10^^ x M'^ erg g-^ M, > 7TeV. (10) 

3. Graviton production through photon fusion and energy loss rate 

Our aim is to study the energy loss mechanism of SN1987A by graviton emission by 
photon-photon annihilation in the ADD framework. For this we need to compute 
the cross-section for the relevant process. Here we present the general formalism for 
calculating the cross-section [19] for two particle initial state. The scattering cross 
section is given by 

1 1 fyj d^pf 



a 



n 



Vrei4:EiE2 J y (27r)32E/ 
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with 



^J{pl ■P2Y - mf 



2^2 



a 



j ^6iE, + E2-E')\M\'. (13) 



Vrel = ^^"^ 7^ (12) 
-C/l-C/2 

where Pi and Ei being the 3-momenta and the energies of the initial particles whose 
masses are mi and Pf and Ej are the 3-momenta and the energies of the final 
particles and /j is the Feynman amplitude for the process. 

For a general reaction of the kind a + 6 — > c, in the center-of-mass frame, the 
expression ffTTj) takes the form 

1 f d^p' 

We use the center-of-mass frame, where we use the following notions. 

= El + E2, (14) 

EiE2Vrel = pVs, (15) 

where p = Pi + P2- 

Next, we focus on the energy loss due to KK gravitons escaping into the extra 
dimensions. The energy loss as per unit time per unit mass of SN in terms of the 
cross-section (7a+b^c, is given by [20] 

{nanhau+b^c)VreiEc) 

ea+b^c. = (16) 

P 

where the brackets indicate thermal average, ria^b are the number densities for a, b and 
p is the mass density and Ec is the energy of the particle c. 

We calculate the cross section using the helicity method [21]- [36]. We follow the 
conventions and Feynman rules derived in [37]. In the helicity method, it is more 
convenient to work with polarizations explicitly. Thus, the polarization vectors [SS] of 
a massive graviton are 

g±2 ^ 2e^e^ 

ejj = v^(eje° + eje^) , 

e% = \j\ + e-et - 26je°) . 

Here and e° are the polarization vectors of a massive gauge boson; for a massive 
vector boson with momentum p^ = [E, 0, 0,p) and mass m, 

e+(p) = -^(0,1,2,0), (17) 
6;(p) = i=(0,-l,z,0) , (18) 
4{p) =-{p,0,0,-E) . (19) 
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The graviton polarization vectors satisfy the normahzation and polarization sum 
conditions 

e^/^'^e^'; = , (20) 

^ ] ^fiu^pa ~ Bfiupa , (21) 

s 

where B^vpa is given by 



+ 2 ( - —r \^^p-^ 



The total squared amplitude, averaged over the initial polarizations z and summed over 
final states for the reaction a'^{qi) + b^' {q2) — * c^" is given by 

(a^(gi) + h^\q2) - c^"(p)) f (23) 

h,h',h" 

where h,h',h" are the helicities and qi,q2,P are the momenta of particles a, b, c 
respectively. 

Photons are quite abundant in supernovae. Here we consider photon-photon 
annihilation to KK graviton and the process is given by, 

l{ki)+^{k2)^KK{p). (24) 

The vertex function for the process is given by [37] 

i 



X 



2Ma 



^al3klfj.k2u — Vp.aklf3k2u — ^ul3klfj,k2a 



+ 1^^^a'l^uf3{kl ■ k2) - ^r]^^ {rjapikl ■ ^2) - kif3k2a) 

+ mnmn-mi-n^iaflvp " ^VfiuVc^f^) + (« ^ /?) • (25) 

The momentum vectors for this reaction are 

p'^ = (m„,0,0,p) (26) 

k'; = {kuO,0,ki) (27) 

k^ = {k2,0,0,k2). (28) 

In helicity formalism the reaction ( IMl) can happen in two ways 

j^{k,)+^^{k2)-^KK^\p) (29) 

7±(A;i)+7^(A;2)--KirO(j9). (30) 

Next, consider these two reactions separately and find their corresponding amplitudes. 
For the reaction described in ( !29l) . the helicity amplitude for the KK graviton emission 
by photon-photon annihilation is 

\M (7^(g) + iHq) - KK^\p)) I = X^^-Pe^{kmq)e^l*{p). (31) 
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1 

"2' 



.±2 
-11 



,±2 _ „±2 _ 1 
-12 "^^21 ~ 2 ' 



(32) 
(33) 



4' = -\- (34) 
The non-zero components of the vertex function are 

vim vl212 vl221 ^2112 v2121 ^2222 ^2211 ^1122 
,vY ,VY ,vY , — , — 

Each of them equal to 



-IK 



h ■ k2. (35) 



2 

where ki ■ k2 = 'm'i/2. 

Substituting the various quantities that we have calculated above in equation ( 13T]) . 
we get 

1^ {l^{q) + iHq) - KK^\p)) I = ^ (36) 
The helicity amplitude for the reaction (1301) is, 

1^ [iHq) + 7^(g) - KK\p)) I = X'^^"^e±(fc)ej(g)ej:(p). (37) 
The polarization tensors for gravitons are given by 

e?i = e°2 = -y| (38) 
e?2 = e^i = 0. (39) 

The non-zero components of the vertex function are X^^^^ — X^^^^ and — X^^^^ 

and are equal to ( l35l) . 

Substituting the various quantities that we have calculated above in equation (!37|) . 
we get 

l-M (7^(g) + l^{q) ^ ^i^°(p)) I = 0. (40) 

Thus the total squared amplitude, averaged over the initial three polarizations and 
summed over final states, is 

^ E 1-^ {l\k,)+i'\k,) KK^'\p))^ = (41) 

h,h',h" 

Substituting this in (fT3|) and using (fT^ and (fT5|) . the cross-section for the process is 
obtained as 

Tlfi^y/S 



a 



-5imH - v^) , (42) 



16 

where s is the center of mass energy, and the mass of the KK state at level n. 
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Since for large R the KK gravitons are very light, they may be copiously produced 
in high energy processes. For real emission of the KK gravitons from a SM field, the 
total cross-section can be written as 

o"tot = K^^a(ra) , (43) 

n 

where the dependence on the gravitational coupling is factored out. The mass separation 
of adjacent KK states, 0{1/R), is usually much smaller than typical energies in a 
physical process, therefore we can approximate the summation by an integration which 
can be performed using KK state density function [37] . 

= (4.)-/T(,,/2) - ("^l 
The volume emissivity of a supernova with a temperature T through the process under 
consideration is obtained by thermal-averaging over the Bose-Einstein distribution 

, 2d% 1 f2d% 



X 



(2vr)3 e^i/^ - 1 J (27r)3 e^^/r _ i 
— y a^^^kk{s,mn), (45) 

n 

where the summation is over all KK states, and the squared center of mass energy s is 
related to the photon energies io\ and 102 and the angle between the two photon momenta 
d^^ as follows: 

s = 2c<JiCt;2(l — cos^^^^) . (46) 
After carrying out the integrals and the summation over KK states, we find 



2'^+3r(f + 3)r (f + 4)C(f + 3)C(f + 4) T"+^ 
(n + 4)7r2 



where we have used M^^'^R"'Sn = M^i and numerically, these Riemann zeta-functions 
are close to 1. In this calculation, we have neglected the plasma effect, through which 
the photons can have different energy dispersion relations from those of free particles. 

We take the supernova core density ~ 10^^ g cm~^. Using ( |T6|) . we compute the 
energy loss rate for n = 2 and n = 3 extra spatial dimensions and hence the lower limits 
on using the conservative upper limits on the energy-loss rate of SN [1] 

e^^ ~ 10^^ erg g-^sec'^ (48) 

The results are summarized below, 

n = 2, e = 4.7 X 10^^ x M"^ erg g'^ sec-\ M, > 14.72TeV, (49) 

n = 3, e = 1.1 X 10^° x M"^ erg g'^ sec~\ M, > 1.62TeV. (50) 

We now combine the energy loss rate due to photon fusion process with that of 
the nucleon-nucleon brehmstralung and rederive the constraints as follows. Cullen and 
Perelstein [5] 

n = 2, e = 6.837 x 10^^ x M"^ erg g-^ s-\ M, > 50.13TeV; (51) 
n = 3, e = 1.131 X 10^^ x M"^ erg g-^ M, > 4.08TeV. (52) 
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Barger, Han, Kao and Zhang [H] 

n = 2, e = 6.747 x 10^^ x M"^ erg g'^ s~\ > 50.96TeV; (53) 

n = 3, e = 6.410 x 10^^ x M"^ erg g-^ M, > 3.64TeV. (54) 

Hanhart et. al. |8j 

n = 2, e = 9.710 x 10^^ x M"^ erg g-^ s'^, M, > 31.39TeV; (55) 
n = 3, e= 1.680 x 10^^ x M"^ erg g'^ M, > 2.79TeV. (56) 

Hannestad and Raffelt P [ID] 

n = 2, e = 4.985 x 10^^ x M"^ er^ M, > 84.03TeV; (57) 

n = 3, e = 1.681 x 10^^ x M"^ erg g'^ s'\ M, > 7.00TeV. (58) 

As we expected, the energy loss rate due to nucleon-nucleon brehmstralung is 1 to 3 
orders of magnitude more than that due to photon fusion process. Hence the combined 
bounds add only the second decimal place to M*. 

4. Conclusions 

We have revisited the constraints on ADD model coming from photon annihilation into 
KK graviton in SN cores. For a conservative choice of the core parameters, we obtain 
the two photon process bounds on the fundamental Planck scale M^, > 1.6 TeV. The 
energy loss rate due to nucleon-nucleon brehmstralung is 1 to 3 orders of magnitude 
more than that due to photon fusion process. Hence the combined bounds add only 
the second decimal place to M^. Thus the present study can strengthen the results 
which are available in the current literature for the graviton emission from SN1987A. 
Our results show that the above processes put a very strong constraints on models with 
large extra dimensions for the case of n = 3. Notice that the plasmon effects are not 
considered in our calculations and will be done elsewhere. 
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